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Intrinsic noise is known to be ubiquitous in Josephson junctions. We investigate a voltage biased 
superconducting tunnel junction including a very small number of pinholes - transport channels 
possessing a transmission coefficient close to unity. Although few of these pinholes contribute very 
little to the conductance, they can dominate current fluctuations in the low-voltage regime. We 
show that even fully transparent transport channels between superconductors contribute to shot 
noise due to the uncertainty in the number of Andreev cycles. We discuss shot noise enhancement by 
Multiple Andreev Reflection in such a junction and investigate whether pinholes might contribute 
as a microscopic mechanism of two-level current fluctuators. We discuss the connection of these 
results to the junction resonators observed in Josephson phase qubits. 

PACS numbers: 73.23.-b, 74.40.+k, 85.25.Cp 



I. INTRODUCTION 

Implementing qubits using superconducting circuits 
[TJ [2] seems to be one of the most promising approaches 
to design a quantum computer. Various implementation 
schemes have been developed, [3 HI El El El [3 El HD] The 
crucial and indispensable device in all these setups is a 
Josephson tunnel junction. Hence, microscopic under- 
standing of this kind of junction and all possible details 
is essential to advance this field. Current fluctuations in 
Josephson Junctions, as they are discussed in this pa- 
per, are of particular importance as they contribute to 
decoherence. 



A. Decoherence, 1/f noise 

One of the major challenges for the realization of prac- 
tical quantum computing is to perform a sufficient num- 
ber of quantum manipulations within the coherence time. 
The need to maintain quantum coherence during the op- 
eration is especially difficult to achieve in solid state sys- 
tems which couple relatively strongly to uncontrollable 
environmental degrees of freedom, that generate quick 
decoherence. 

After electromagnetic qubit environments have been 
successfully engineered to improve coherence, we are now 
mostly concerned with intrinsic noise of the solid state 
system. The most prominent source of intrinsic decoher- 
ence is non-Gaussian 1/f noise [TT], for which the spectral 
function behaves like S(u>) oc 1/uj [T2][T3]. 1/f noise typi- 
cally appears due to slowly moving defects in strongly dis- 
ordered materials and is usually explained by an ensem- 
ble of two-level fluctuators coupling to the system under 
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consideration. A heat bath causes uncorrelated switching 
events between the two states, which are described by a 
Poissonian distribution with mean switching time r. For 
a single fluctuator this leads to random telegraph noise 
(RTN). Superimposing several such fluctuators, using an 
appropriate mean switching time distribution p(r), re- 
sults in a 1/f noise spectrum. 1 // noise seriously limits 
the operation of superconducting qubits [TH [TS] . 

Besides magnetic-flux fluctuations [T5J [17], critical- 
current fluctuations due to charge trapping at defects 
in the tunnel barrier [TH] is a prominent, possible mech- 
anism for low- frequency 1/f noise in junctions of super- 
conducting qubits. As complement to this, we will inves- 
tigate the intrinsic noise of Josephson tunnel junctions 
containing a few high-transmission channels, that poten- 
tially reside in the junction. We will address the question 
whether such defects might introduce another intrinsic 
source of current fluctuations leading to 1/f noise. 



B. Rough superconducting tunnel junctions 

A Josephson tunnel junction consists of two supercon- 
ductors separated by a thin insulating oxide layer. Trans- 
port through such small scatterers can be described by 
quantum transport channels [19 which, in our case, refer 
to the channels of the oxide. 

In tunnel junctions, the transmission of all transport 
channels are assumed to be small. However, the fab- 
rication process is usually not at all epitaxial quasi- 
equilibrium growth, thus one has to expect the oxide 
layer to be non-crystalline and disordered [201 HI]- We 
will take this into account by investigating rough super- 
conducting tunnel junctions, where we assume that the 
junction additionally contains a few transport channels 
with very high transmission eigenvalues. Pinholes, see 
Fig. [l] might occur as defects due to the fabrication pro- 
cess. Indeed, in Ref. [22], the importance of pinholes was 
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FIG. 1: Schematic diagram of the oxide layer of a Joseph- 
son Junction. Several transport channels are indicated. The 
arrow thickness corresponds to the respective transmission 
eigenvalue. Left: schematic diagram of an ideal Josephson 
tunnel junction. The transmission of all channels is small. 
Right: considered rough Josephson tunnel junction, i.e., a 
tunnel junction including some high-transmission channels, 
so-called pinholes. 

pointed out, but also Ref. [S3] discusses the significance 
of pinholes in mesoscopic devices, e.g., for the Kondo ef- 
fect. 

There has been particular interest in pinholes to un- 
derstand subharmonic gap structure in weak links. In 
Ref. |24j the subharmonic gap structure of a tunnel junc- 
tion was modeled by assuming that some channels have 
pinhole character. Analyzing superconducting qubits 
containing pinholes has additional motivations that are 
going to be reviewed in the following subsections. 

1. Multiple Andreev Reflections (MAR) 

In Josephson Junctions with voltage bias smaller than 
the superconducting gap, direct tunneling of quasipar- 
ticles is impossible. In this case, charge transport is 
governed by Multiple Andreev Reflection (MAR), see 
Ref. [251 [261 [27] . 

Andreev Reflections occur at contacts between a scat- 
terer and a superconductor. In a system consisting of 
two superconductors separated by a scatterer there are 
two superconductor-scatterer interfaces where Andreev 
Reflection can take place. This leads to processes in- 
volving a sequence of Andreev cycles known as Multiple 
Andreev Reflection (MAR), in which charge can be trans- 
ferred even for voltages smaller than 2A/e. In general, 
an n th order MAR process, transferring n charge quanta 
at a time, a so-called Andreev cluster, comprises (n — 1) 
Andreev Reflections and occurs above a threshold volt- 
age V n — 2 A / (en) . For voltages below this MAR voltage 
V n , the energy gap cannot be overcome by (n — 1) re- 
flections. As these processes are composed of multiple 
transmission cycles, they sensitively depend on the elec- 
tron transmission probability, i.e., the set of transmission 
eigenvalues characterizing the junction. 

We can thus expect that rough superconducting tunnel 
junctions will be highly affected by MAR and we will see 
that even very few pinholes will have an extreme impact 



on the transport properties of the junction. 

2. Noise enhancement due to MAR 

In Ref. 22], shot noise 28] of NbN/MgO/NbN 
superconductor-insulator-superconductor tunnel junc- 
tions was measured. The result shows noise enhanced 
by Andreev reflection. The authors attribute this to the 
occurrence of MAR processes in pinholes in the MgO bar- 
rier. They model their data assuming Poissonian shot 
noise 2el, where they replaced the single charge quan- 
tum e by an effective transferred charge q(V), due to 
MAR. 

Such processes might be highly relevant as a source of 
intrinsic noise in superconducting qubit devices due to 
pinholes residing in the Josephson junction. It has to be 
realized that, in the case of transport through pinholes, 
the shot noise is not governed by the simple Poisson for- 
mula which is only valid in the limit of small T. The 
method, we will use, properly treats all possible trans- 
mission eigenvalues, and quantitative statements on the 
impact of rough barriers will be possible. 

3. Junction Resonators 

A new measurement revealing major intrinsic sources 
of decoherence in Josephson junction qubits was per- 
formed in Ref. [21] . The authors observed characteristics 
of energy-level repulsion at certain frequencies, as pre- 
dicted for coupled two-state systems. This structure of 
level-splittings was attributed to spurious resonators re- 
siding in the Josephson junction. Measurements of Rabi 
oscillations revealed that these resonators cause signifi- 
cant decoherence [3D]. Similar to the scenario of charge 
trapping, mentioned before with respect to 1/f noise, the 
energy-level repulsion could be explained by assuming 
two-state current fluctuators in the junction. 

Although other processes such as charge trapping 
within the junction barrier are believed to be relevant 
effects for realizing such spurious resonators, pinholes in 
rough tunnel junctions might be additional candidates 
for introducing two-state current fluctuators, see Sec- 
tion EE 

The structure of this paper is as follows: After a short 
survey concerning the method used, we discuss leakage 
current of rough superconducting tunnel junctions. This 
is followed by a section regarding its noise properties. 
Finally we discuss the full counting statistics of pin- 
holes and whether they might contribute as a microscopic 
mechanism of two- level current fluctuators. 



II. METHOD 

We will be interested in the full counting statistics 
(FCS) of charge transfer through the junction, i.e., the 
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probability distribution P to (N) for N charge quanta to 
be transmitted within measurement time to [311 132[ 133] . 
In addition to the noise characteristic, proportional to 
the second cumulant, this distribution will also supply 
us with higher, non-Gaussian cumulants such as they oc- 
cur in RTN. 

To calculate P to (N) we will apply the non-equilibrium 
Keldysh Green's function approach 34J. Within this 
scheme, it is possible to employ several quantum field- 
theoretical methods, used in the transport theory of 
metals |35j . and describe the system under considera- 
tion microscopically. In the zero temperature limit, the 
cumulant generating function of a single mode voltage- 
biased Josephson Junction was calculated analytically in 
Ref. [361 157] . We extend this result to multimode junc- 
tions containing M transport channels which are char- 
acterized by a set of transmission eigenvalues {T n }. In 
this case, based on [36l[37], we immediately find for the 
cumulant generating function 



eV 



dE 



x In 



l + J2 p n(E,V,T n )(e m X-l) 



n=0 



, (1) 



where P n (E,V 7 T) is computed according to Ref. [57] . 

With this approach, we accurately take into account 
MAR within the rough junction model specified above. 
Thus we will find quantitative results going beyond an 
effective charge q(V) alone. 



III. LEAKAGE CURRENT 

We quantitatively investigate leakage current, i.e., cur- 
rent in the subgap voltage regime eV < 2A of a voltage- 
biased rough superconducting tunnel junction. Using 
Eq. ([lj, the average current of our junction, containing M 
transport channels characterized by a set of transmission 
eigenvalues {T„} which is described by the distribution 
p(T), is given by 



I = 



2e 



M \ dTp(T) I dEJ2nP n (E,V,T). 
Jo Jo 



(2) 



A. Homogeneous multimode contacts 

For illustrative reasons we start from a homogeneous 
multimode contact between superconductors containing 
M transport channels all with the same transmission 
eigenvalue T\. The transmission eigenvalue distribution 
reads p(T) = S(T - T x ). In the case M = 1, this would 
be a single-mode quantum point contact (QPC). We com- 
pute conductance in units of the normal state conduc- 
tance G N = J dTp(T). From Eq. (B it is clear 
that the normalized average current of a homogeneous 




1e-15 - 



1e-20 



0.5 1 1.5 2 2.5 

eV/A 



0.5 



1 1.5 

eV/A 



2.5 



FIG. 2: Current through a homogeneous multimode contact 
with small transmission between superconductors as a func- 
tion of bias voltage at T — OK on a logarithmic current scale. 
Inset: linear current scale. 



multimode contact is the same as the one of a single- 
mode QPC which was already discussed in Ref. [37] . 

For small transmission eigenvalues Fig. [2] shows the av- 
erage current as a function of bias voltage for low trans- 
mission probability Ti < 1 on a logarithmic and linear 
scale. We see that a contact with T% = 0.1 already devel- 
ops a relatively large leakage current in the subgap regime 
eV < 2A. Another immediate aspect, which will become 
important below, is its scaling. Conductance steps of 
size T\ arise at MAR voltages 2A/n, demonstrating that 
the current is reduced by a factor of T\ at each step. 
For a single-mode QPC it was shown before that cur- 
rent transport for small transmission eigenvalues in the 
voltage interval [2A/(ne),2A/(n — l)e] is dominated by 
the n tft, -order MAR process. In Ref. [37 the authors ex- 
plicitly showed this T n dependence within a perturbative 
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FIG. 3: Current through a homogeneous multimode contact 
with high transmission between superconductors as a func- 
tion of bias voltage at T — OK on a logarithmic current scale. 
Inset: linear current scale. 



4 



calculation. 

For high transmission contacts we note that pertur- 
bative approaches in T will fail and it is necessary to 
use non-perturbative methods, as we do here. This will 
become in particular important for deriving quantita- 
tive results for rough junctions containing low- and high- 
transmission channels. Fig. [3] shows the current for a 
range of transmission probabilities T > 0.6. Especially 
at small voltages, the current through high-transmission 
modes is larger by orders of magnitude compared to the 
small transmission case. 



B. Rough tunnel juncitons 

We now turn to rough Josephson tunnel junctions as- 
suming a small number of pinholes with transmission 
eigenvalues close to unity that reside in the junction. We 
consider a contact with M channels. A fraction a of these 
channels has a high transmission eigenvalue T\ , the vast 
majority has a small value T 2l typical for tunnel contacts. 
Altogether, we consider the eigenvalue distribution 

p(T) = aS(T-T 1 ) + (l-a)6(T-T 2 ), (3) 

causing a normal conductance of 



small pinhole fraction a is shown in Fig. [4] We see 
that, starting at high voltages, the current follows the 
well-known subharmonic gap structure curve for tunnel 
transmission coefficient T 2 only above a certain voltage, 
depending on a. However, below that voltage, due to 
the fact that the current carried by tunnel transmis- 
sion channels is reduced by a factor of T 2 each time the 
voltage passes another MAR voltage 2A/n, the highly 
transmissive channels dominate, leading to a smooth, 
weakly structured subgap contribution that level off into 
a plateau before it drops again. Consequently, at suf- 
ficiently low voltages, the current through rough tunnel 
contacts is overwhelmingly carried by the pinhole frac- 
tion. 

With respect to the pinhole fraction a, our quantitative 
treatment demonstrates that in tunnel junctions we can 
only have very few pinholes, which Refs. |24j seems to 
overestimate. From Fig. |4j for example, we infer that 
in a junction, where a possible measurement shows two 
full current steps each scaling with a factor of T = 0.01 
at eV — 2 A and eV — 2 A/2, respectively, we can have 
roughly less than 1 out of 10 6 channels with transmission 
T > 0.6! 



C. Characterizing Pinhole Thresholds 



Gn — 



r A/[aTi + (l-o)T 2 ]. 



In Section [ill C| we are going to discuss why this distribu- 
tion captures the essential physics of more complicated 
distributions. 




FIG. 4: Leakage current on a logarithmic scale as a func- 
tion of bias voltage at T = OK for a rough superconduct- 
ing tunnel junction with transmission eigenvalue distribution 
p(T) = a5(T-T 1 ) + (l-a)S(T-T2), Ti = 0.6, and T 2 = 0.01. 
Inset: Ti = 0.6, and T2 = 0.003. The different curves refer to 
different pinhole fractions a within all transport channels. 

We calculate the current for this model, taking into 
account the very different transport properties of T\ and 
T2 channels. The result for two scenarios with extremely 



As we have seen, below a certain voltage, a high trans- 
mission channel residing in a rough Josephson tunnel 
junction dominates high-order subharmonic steps in the 
current characteristic. We can use this result to charac- 
terize the fraction of pinholes in all transmission channels 
by very sensitive current measurements. 

In Ref. [SB] a current- voltage plot for an Al-Al 2 0^-Al 
junction used in a Josephson-j unction qubit is presented. 
At eV = 2A the measured current shows a subharmonic 
step corresponding to a tunnel transmission eigenvalue 
of T = 0.003. The second current drop at eV = A is 
indicated but the measurement does not resolve the next 
expected plateau. The experimental result is consistent 
with the calculation presented in the inset of Fig. [4] for 
a pinhole fraction a = 0.5 • 10~ 5 . This corresponds to 
one pinhole of Ti = 0.6 in a junction of 1/a = 200 000 
channels. Actually, in Ref [29], the number of transport 
channels for the junction under consideration [38] is es- 
timated to this order of magnitude, indicating that the 
existence of pinholes in state-of-the-art superconducting 
qubit devices is compatible with current measurements. 

Indeed, new design concepts have lead to a significant 
reduction of the junction size (see Ref. [39 ) and with 
it a supression of intrinsic noise. Clearer insight would 
be provided by highly sensitive current-voltage measure- 
ments at voltages stretching out over several current 
steps at V n = 2A/n. In the following, we will assume 
a small number of pinholes. 

Finally, this provides the justification for the very sim- 
ple transmission distribution function ^ . The transmis- 
sion eigenvalues are determined by WKB, T = exp(— nd), 
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and this way, depend on the junction width d. Then 
the pinhole transmission eigenvalues might be related 
to a distribution if widths p(d) of the oxide layer sep- 
arating the superconductors. Considering the strict non- 
negativity of d, a lognormal distribution might be appro- 
priate for describing p(d) for the pinholes. 

All this can be done in our approach, but as we have 
seen above, in state-of-the-art superconducting qubit de- 
vices we might only have a small, single-digit number of 
pinholes in a huge junction. Thus, doing statistics is not 
necessary and considering a single value Ti to represent 
the pinhole transmission eigenvalue distribution, as done 
in Eq. ([3]), is a sufficient way to take them into account. 



IV. NOISE 

We will examine the noise properties of rough super- 
conducting tunnel junctions. As we have seen a small 
amount of pinholes contributes very little to the conduc- 
tance in the low-voltage regime. Now we are going to 
show, that pinholes on the contrary do dominate current 
fluctuations. 

Using the cumulant generating function (JlJ, we find 
for the zero-frequency noise 

•'° n n 

(4) 

A. Homogenous multimode contacts 



el / G N A 
1 S,/(2G N A) 



Again, as in Section we start from homogeneous 
multimode contacts with each channel having the same 
transmission eigenvalue and begin with small transmis- 
sion T< 1. Fig. [5] shows noise and current characteris- 
tics. In this there is always one very dominant 
MAR- process, which causes charge transport to be dom- 
inated by quanta of q(V) = e(l + [2A/eUJ), the noise 
scales with this additional charge factor. Thus, in the 
small-transmission regime Poissonian shot noise Si = 2el 
with modified charge quantum nq(V), properly explains 
the observed noise features. 

In the case of large-transmission eigenvalues, inset 
of Fig. [5j the noise characteristic changes dramatically. 
For very high probabilities T, the noise increases with 
decreasing voltage in the subgap regime. Depending 
on the value of T, it develops a maximum, but falls 
off again at even lower V. Remarkably, and in strong 
contrast to any simpler model, we note that a con- 
tact with perfect transmission T = 1 shows low, but 
finite noise. This is markedly different from the nor- 
mal conducting case where, given the shot noise formula 
Sj = (Ve 3 /irh)T(\ — T), we would anticipate zero noise 
in the case of perfect transmission. Furthermore, we see 
that the larger the transmission the steeper and higher 
is the noise ascent for small voltages. For high eV, the 
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FIG. 5: Noise and current as a function of bias voltage for 
a homogeneous superconducting tunnel junctions with trans- 
mission eigenvalue T = 0.01. Inset: noise characteristics for 
junctions with large transmission. 
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FIG. 6: Noise Si as a function of transmission eigenvalue T 
for a homogeneous superconducting tunnel junction where the 
bias voltage is set as parameter. Inset: enlargement around 
T = 1. 



high transmission curves approach the T = 1 charac- 
teristic. Thus, altogether we see that in this case the 
description with pure Poissonian shot noise with modi- 
fied charge quantum is insufficient and the generalization 
used in the rest of this paper shows new features. 

It is instructive to look at the noise curve from a dif- 
ferent perspective. Focussing on the T-dependence, in 
Fig. [6j we set voltage as a parameter and plot noise as 
a function of transmission. The noise develops a maxi- 
mum at high transmission values. As noticed before, each 
curve falls off to a finite residual noise level at T = 1. 
For smaller voltages, the maximum becomes more and 
more pronounced and it seems to be squeezed into the 
high-transmission regime. On the order of eV = 0.1 only 
channels with very high transmission significantly con- 
tribute to the noise. 

The explanation of the noise features for high transmis- 
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sion is more involved. From Eq. Q we see that the noise 
can be expressed in terms of the variance of P„(P, V, T), 
which is the probability for a MAR-process transferring 
charge quanta ne 



Si 



4c 2 

IT 



dE[(n 2 



4c 2 

IT 



dE Var(n) 



( 5 ) 

Thinking of shot noise as partition noise, for a single- 
mode normal conductor with perfect transmission T = 1, 
there is no uncertainty whether a particle is transmitted 
or reflected. We find zero noise. In the superconduct- 
ing case, due to perfect transmission, we are still certain 
about charge transfer taking place, but an additional un- 
certainty is introduced. For high transmission includ- 
ing T = 1, there are many different MAR-processes con- 
tributing to charge transport, which is described by the 
probability distribution P n . This additional uncertainty 
is the qualitative physical explanation of the finite noise 
observed in the case of perfect transmission. 



1. Toy Model 

To clarify the essential physics, referring to the full 
computer-algebraic calculation is unsatisfactory. Thus, 
we will try to explain the basic noise features with the use 
of the toy model presented before in Ref. [37]. Originally, 
this model was introduced to illustrate how to calculate 
the cumulant generating function of a weak link with 
voltage bias in an easy, analytically solvable case. We 
summarize the basic simplifying assumptions. 

We only look at voltages equal to one of the MAR 
voltages eV = 2A/n, and for each of them we only take 
into account one MAR process, namely the most relevant 
one which transfers 



N = 



2A 

eV 



1 



(6) 



charge quanta. This simplifies the cumulant generating 
function S(x) to the one of a binomial distribution. Fur- 
thermore, in this model, Andreev reflection above the 
gap is neglected and the Green's function is simplified by 
assuming a constant density of states above the gap. 

The cumulant generating function for the toy model 
reads 



S(x) 



2eVt 



ln[l + P w (e w *_i)] 



(7) 



For every MAR voltage V, another specific transport 
process with probability Pn is relevant, Eq. We 
emphasize that, due to this, the argument of the loga- 
rithm in Eq. ^ depends on voltage via the selection of 
the relevant Pn- Fig. [7] shows the toy-model probabili- 
ties Pjv as a function of transmission eigenvalue T. For 
perfect transmission, as we have reduced the system to 
a binomial distribution involving only a single transport 
process, each probability is unity. The probabilities for 



N > 2 and imperfect transmission are always smaller 
than in the normal conducting case, because a higher- 
order process is necessary in order to transfer charge. 
For large charge quanta N, very high transmission is 
necessary, since many Andreev reflections are involved 
in such a process. Thus, for larger N, i.e., small voltage 
bias (Eq. ([6])), nonzero probabilities are more and more 
shifted to the high-transmission regime. 
The second cumulant 



C 2 = N- 



2eVt ( 



Pv(l-Pv) 



(8) 



is proportional to the noise correlator. The expression 
Piv(l — Pn), which matches the one in the traditional 
shot-noise formula if we replace T by Pn, is displayed 
in the inset of Fig. [7] For large N, i.e., small volt- 
age, the maximum is shifted and squeezed into the high- 
transmission regime. 

So, altogether, we distinguish two mathematical ingre- 
dients to the noise. One is the expression Pjv(1 — Pzv) 
that we just discussed. Additionally, there is the prefac- 
tor N 2 (2eVt /h). In the small voltage regime it results 
in noise enhancement that behaves approximately like 
1 /V. As the noise is determined by the product of both 
parts, for a fixed transmission coefficient, there will be a 
voltage regime where the noise gets enhanced by lowering 
the applied voltage bias. However, at some voltage, or 
conversely for some N, Pjv(1 — P/v) will overcompensate 
this increase and reduce the noise again. To summa- 
rize, the toy model still explains noise enhancement by 
an increased charge quantum. The decrease of noise at 
very low voltage follows from the overcompensation of 
this effect by the decrease of transfer probability in the 
expression P/v(l — Pat). 

In the inset of Fig. [8j for comparison, the noise calcu- 
lated using the full expression and the toy model at MAR 
voltages, is presented in a single plot. The simplified 
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FIG. 7: Toy model probability Pn to transfer N charges at 
a time as a function of transmission eigenvalue T. Inset: ex- 
pression Pjv (1 — -Pjv) that occurs as a term in the 2 nd cumulant 
(Eq. |8]|) as a function of transmission eigenvalue T. 
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FIG. 8: Maximal noise maxT [Si(T, V)], optimized with 
transmission as a parameter, as a function of bias voltage on a 
double logarithmic scale for a homogeneous contact between 
superconductors. The calculated data using the full expres- 
sion, Eq. 0, are fitted using the given power law. Inset: 
Comparison of the noise results calculated for a T — 0.97 ho- 
mogeneous contact between superconductors as a function of 
bias voltage using either the full expression or the toy model. 



model qualitatively shows the basic features of our nu- 
merical calculation. Nevertheless, there is a huge quanti- 
tative difference. Thus, we realize that the toy model is 
qualitatively sufficient but it fails dramatically to provide 
quantitative results. Thus, for quantitative calculations, 
the extensive calculation used in the remainder of this 
paper is essential. 



2. What is driving the noise increase? 



Noise of rough superconducting tunnel 
junctions 



We now return to the model of Section |fff C| There, 
we looked at a rough superconducting tunnel junction 
with eigenvalue distribution given by Eq. ([3]). Here, we 
are concerned with the noise generated in this kind of 
junction. Fig. [9] shows the result. 

In contrast to a normal tunnel junction we see a dra- 
matic change in the noise characteristic due to very few 
pinholes with an enormous noise increase at small volt- 
ages. As in the case of leakage current, at a certain point 
in the subgap regime, the pinholes begin to dominate the 
noise characteristic. In this range, Si is solely carried by 
the few pinholes. 



Together with our results in Section |III| this demon- 
strates one of our central results: although a small 
amount of pinholes residing in the junction contributes 
very little conductance, it can dominate current fluctua- 
tions in the low-voltage regime. As pointed out before, 
sensitive measurements of the leakage current will pro- 
vide an estimate on the amount of pinholes that might 
be contained in the considered junction. 

The considered pinhole transmission eigenvalue of T = 
0.986 is chosen in order to display all structure at volt- 
ages down to eV = 0.05 A. Nevertheless, analogous to 



Section [IV A| we can add two more aspects: Firstly, for 
smaller voltages than resolved in Fig. [9] the noise will 
show a maximum and then will fall off again. Secondly, 
considering higher values of transmission will lead to an 
even steeper and higher ascent, starting at smaller volt- 
ages. 



We can ask the question: what is the maximal noise 
at a given voltage? This means, for fixed voltage bias, 
we use the transmission eigenvalue as a parameter to 
find the maximal value. In the toy model, Fig. [7] 
maxy [-Pjv(l — Pn)] is always 1/4. Thus, here, the maxi- 
mal noise maxy [Si(T, V)] depends only on the prefactor 
in Eq. ([8]). Consequently, for small voltages, it approxi- 
mately scales like l/V. 

For the full theory, in Fig. [8] the maximal noise 
maxT [Si(T, V)] is plotted against voltage bias on a dou- 
ble logarithmic scale. In the small-voltage regime, the 
data can be fitted well using a power-law. We find 

max[Si(T, V)]<x^j. 

Thus, although quantitative statements resulting from 
the toy model and from the full expression differ signif- 
icantly, we see that the maximal noise at given voltage 
follows a similar power law with an exponent of 0.8 in- 
stead of unity. Hence, even in the much more compli- 
cated situation, including multiple MAR processes, the 
inherent l/V dependence, which basically results from 
the increased charge quanta, seems to play a major role. 
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FIG. 9: Noise Si on a logarithmic scale as a function of bias 
voltage for a rough superconducting tunnel junction charac- 
terized by the transmission eigenvalue distribution p(T) = 
a6{T-T 1 ) + (l-a)8(T-T 2 ), Ti = 0.986 and T 2 = 0.01. The 
different curves refer to different pinhole fractions a within all 
transport channels. Inset: linear noise scale. 
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FIG. 10: FCS for a transport channel between superconduc- 
tors with transmission T = 0.936 and voltage bias eV = 0.3A 
at measurement time to — 10 Tj. Inset: measurement time 
to = 100 Tj 



V. FULL COUNTING STATISTICS OF 
PINHOLES 

We will investigate the full counting statistics (FCS) 
of charge transport through pinholes of a rough super- 
conducting tunnel junction. Given the definition of the 
cumulant generating function, the FCS is calculated by 
Fourier transformation. Using this insight into the trans- 
port process, we will furthermore discuss a possible model 
of high transmission channels as microscopic origin of 
two-level current fluctuators. 



A. Resolution of structure in the Full Counting 
Statistics 

To determine the probability distribution P to (iV), we 
have to set the measurement time t . In general, the 
calculation of the cumulant generating function for a 
voltage biased Josephson junction is significantly compli- 
cated due to the ac Josephson effect [36] [37] . In order to 
make computation feasible, and to avoid interpretation 
difficulties of arising 'negative probabilities' in the su- 
perconducting system [40], to must be sufficiently longer 
than the inverse of the Josephson frequency Tj = h/2eV. 
(See Ref. [37] for further details). Consequently Tj sets 
a time scale in our approach and there is a lower bound 
for the measurement time to- 

We consider a contact with transmission eigenvalue 
T = 0.936 at low bias voltage eV — 0.3A, where qubits 
might be operated, and take into account two different 
measurement times t = 10 Tj and t = 100 Tj. Fig. [10] 
shows the results: For the long measurement time the 
FCS is Gaussian. In contrast, for to = 10 Tj we see a 
rich comb structure. 

We will discuss this comb structure and its origin in 
detail later on. Here, we want to point out that this 



structure turns into a Gaussian for long measurement 
time to- This is as we would expect: if we sum the num- 
ber of transfered charges over a very long time it will 
become possible, instead of considering individual MAR 
processes with their specific probabilities, to just assign 
an average likelyhood for one elementary charge quantum 
to be transfered. Thus, in the long measurement time 
limit, transport can be described by a sum of many inde- 
pendent and identically-distributed events what results 
in a Gaussian. This is the essence of the central limit 
theorem used in statistical physics. Indeed, the prob- 
lem above can be related to the quasi- ergodic hypothesis. 
Hence, it is clear that for very long measurement times 
the comb structure, due to individual, discrete transport 
processes, is washed out. 

For significantly higher voltages, like eV = 1.5 A, the 
most relevant transport processes transfer much smaller 
charge quanta. It turns out that consequently, in this 
case, discrete structures in the FCS cannot be resolved 
using the time interval to = 10 Tj. Despite these limita- 
tions concerning to, we can resolve structure in the FCS 
for a limited parameter window. 

B. Pinholes as Junction Resonators? 

1. Motivation 

We are now coming back to the questions whether pin- 
holes might explain decoherence from junction resonators 
in phase qubits or 1/f noise P21 [13j [18]. Thinking of 
the different possible MAR processes, which transfer dif- 
ferent sizes of charge quanta, a pinhole might introduce 
current fluctuators: imagine a high transmission chan- 
nel, i.e., a pinhole hidden in the junction. Two different 
MAR processes A and B transfer charge in two differ- 
ent quanta n^e and rise. Thus, we might think of two 
current states \A) and \B); each of them carry charge 
using only one of the distinct MAR processes A, B, re- 
spectively. Due to the differently sized Andreev clusters 
being transfered, the two states will cause two different 
currents. In principle the mechanism is similar to the idea 
of charge-trapping [18], where a trapped charge blocks 
tunneling through a transport channel. There, one intro- 
duces an untrapped state |r u ) causing high current and 
a trapped state |Tt) causing low current. In comparison, 
we consider two current states \A) and \B) corresponding 
to charge transport by two different MAR processes and 
thereby causing two distinct currents. 

2. Calculation 

We invstigate whether this scenario results from a pin- 
hole model. If this was the case, we would expect to 
find two distinct peaks in the FCS, where the first one 
refers to charge transport due to MAR process A and the 
second one corresponds to MAR process B, each within 
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the time interval to- Hence, let us see whether we find 
parameters that result in such an FCS. 

We consider very high transmission channels, for in- 
stance T = 0.99 and calculate the FCS for this transmis- 
sion eigenvalue at two subgap voltages. The results are 
shown in Fig. We find two very pronounced peaks in 
the FCS. Note that here the measurement time is very 
short but, despite some artifacts in the diagrams, the 
distribution still has a normalization close to unity. 



3. Attempted interpretation in terms of two-level fluctuator 

Given these pronounced peaks, does this result indi- 
cate a scenario where a pinhole via its different MAR pro- 
cesses might actually introduce a two- level current fluctu- 
ator? If we assume so, we associate the first peak with the 
case where charge transport is carried by MAR process 
A, i.e. charge transport in quanta of n^e only. Accord- 
ingly the second peak refers to the case where transport 
takes places via MAR process B, using charge quanta 
n B e. 



Taking a closer look at Fig. 11 reveals a sharp bound- 
ary for the appearance of peaks towards large total charge 
numbers N. In contrast, to the left, i.e. towards smaller 
N, we see small peaks next to the dominating ones. In 
fact, for a given voltage, there is a lower threshold for the 
MAR order, i.e., a lower bound on the minimal charge 
cluster being transferred in a single MAR process. Fur- 
thermore, in addition to the dominant processes A and B, 
there will also be finite probability for MAR of higher or- 
der, i.e., current flow via even larger quanta than n^e or 
n^e. So according to the two-level interpretation, iden- 
tifying each peak with charge transport due to different 
MAR processes, we would expect this boundary to be 
reversed, namely: a sharp boundary for the existence of 
peaks towards small N, due to the lower bound on the 
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FIG. 11: FCS for a transport channel between superconduc- 
tors with transmission eigenvalue T = 0.99 and voltage bias 
eV = 0.290A. Insert: eV = 0.227A. The measurement time 
is t = 4Tj. 



charge cluster size, and additionally, little peaks towards 
large N, due to the finite probability for MAR of higher 
order than the two dominant ones, A and B. 

The second aspect is the spacing between the peaks. 
For 



1 1 the distance 



the distributions in Fig. 
slightly larger than the smallest possible charge quantum 
( L2A/eV^J + 1), i.e., for the main panel 7 (9 for the inset). 
With respect to the MAR threshold, this is roughly the 
size of the average charge quantum that we would expect 
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to be transferred by a single Andreev cluster. In Fig 
from the number of transferred charges and the minimal 
Andreev cluster size, we infer that, within the measure- 
ment time to, roughly 5 MAR processes contributed to 
the rightmost peaks. In the above two- level scenario, A 
and B are adjacent MAR processes meaning their trans- 
ferred charge quanta differ only in one elementary charge. 
Thus, if a pinhole introduced a two- level current fluctu- 
ator where each peak refers to current flow via distinct 
MAR processes A and B, in Figure [TT] we would expect 
a peak spacing of AiV = 5 rather than a value larger 
than (|_2A/eFj +1). This makes the two-level fluctuator 
hypothesis inconsistent. 



4- Alternative, consistent interpretation 

Thus, the structure we have seen in the FCS of a pin- 
hole does not correspond to the scenario of a two-level 
current fluctuator as suggested above. In fact, the de- 
scription of the probability distribution becomes consis- 
tent if we identify each peak with the number of attempts 
being successful to transmit an Andreev cluster: within 
the measurement time to we might think of a total num- 
ber of attempts to transfer charge cluster, where the 
actual size of the quantum might differ due to the in- 
dividual, possible MAR processes. In the distributions 
of Fig. [TTJ each rightmost peak corresponds to the case 
where every attempt is successful to transfer an Andreev 
cluster, so we get the sharp boundary observed for the 
appearance of peaks towards large N. The next peak 
to the left corresponds to the case where exactly one at- 
tempt fails and so on. Thus, the peaks are naturally sep- 
arated by a distance larger than (|_2A/eVj +1), namely 
the average Andreev cluster size transferred in case of 
a successful attempt. As the actual size of successfully 
transmitted clusters might differ due to the individual 
MAR processes, the pronounced peaks are broadened. 
The comb structure in Fig. [TSJcan be explained the same 
way. Here, in contrast to Fig. [TTJ due to smaller trans- 
mission, the case where every attempt is successful is not 
the most likely one. 



5. Conclusion 

To summarize this section, we have discussed the pos- 
sibility of a pinhole to introduce a two-state current 
fluctuator due to its different MAR transport processes. 
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This is conceptually similar to the mechanism of charge- 
trapping, Ref . [18] . Although at first sight it is suggestive 
to relate the observed peak structure to distinct MAR 
processes, a more detailed analysis suggests a very dif- 
ferent but consistent interpretation in terms of successful 
transport attempts of Andreev cluster. Taking this into 
account, we see no clear evidence that a pinhole might 
be a microscopic origin for introducing two-level current 
fiuctuators. Charge-trapping injunctions is probably one 
of the most relevant mechanisms. However, it might be in 
particular interesting to think about such a process open- 
ing and closing a very high transmission channel i.e. a 
pinhole. Due to the large charge quanta being transfered, 
the process of trapping and untrapping might result in 
high magnitudes of current fluctuations. This picture 
may change if electron-electron interaction is included, 
given that the large charge quanta in a pinhole may effi- 
ciently block large parts of the junction. 

A very intuitive picture might be an occupied upper 
Andreev bound state jH] , that causes a repulsion within 
the channel. Nevertheless, in the case of voltage bias, 
such a state with energy Ej = A[l - Tsin 2 (^(f)/2)] 1 / 2 , 
where <p(t) is the superconducting phase, might be adia- 
batically carried above the gap within the actual Joseph- 
son cycle directly after population. Further research 
might clarify this scenario. 



VI. CONCLUSION 

We have investigated voltage-biased rough supercon- 
ducting tunnel junctions containing some high transmis- 
sion channels, pinholes. We have accomplished this us- 
ing the method of full counting statistics formulated 
within the non-equilibrium Keldysh Green's functions 
technique. Based on this microscopic approach, we were 
able to properly quantify physical effects due to low- and 
high- transmission channels in a single junction. 

By exploring leakage current of such systems, we ob- 
served that a tunnel junction may contain much fewer 
pinholes than previously speculated [53]. We further 
demonstrated how highly sensitive current measurements 
can clarify the existence of pinholes. We pointed out that 
existing current measurements done for junctions of the 
superconducting qubit devices [53], do not strictly rule 
out the existence of a hidden pinhole. 

Furthermore, we examined noise properties. We 
demonstrated that even very few pinholes give rise to 
a drastic increase of the noise in the very low subgap 
voltage regime. Thus, although few pinholes might con- 
tribute very little to the average current, they can dom- 
inate current fluctuations. Although details of this noise 
enhancement, comprising contributions of several MAR 
processes, turned out to be quite complicated, we pro- 
posed that the physical essence of the observed noise 
boost still lies in the increased charge quantum that is 
transfered. To do this, we compared the explicit noise 
calculation to a simplified model. This showed qualita- 



tive agreements, and thus illuminated some essential fea- 
tures, but failed quantitatively, therefore demonstrating 
the need of a full calculation. 

Finally, we investigated the full counting statistics 
(FCS) of charge transport through pinholes. Despite lim- 
itations concerning the measurement time to, we could re- 
solve non-Gaussian structure in the FCS for a limited pa- 
rameter window. We discussed a possible model of high- 
transmission channels as a microscopic origin of two-level 
current fiuctuators. Indeed, for certain voltage parame- 
ters, the FCS shows a two-level peak structure. From 
a more detailed analysis we inferred that this structure 
cannot be related to charge transport by distinct MAR 
processes. Thus, given the dc part of the probability dis- 
tribution, we find no evidence that a pinhole might intro- 
duce an additional source of two- level current fiuctuators. 
We presented an alternative, consistent interpretation of 
the observed peak structure in terms of successful trans- 
mission attempts of Andreev clusters. 

So far, our approach is limited to the stationary or 
quasi-stationary case. Improvements on this might incor- 
porate time-dependence into the Keldysh Green's func- 
tion approach. This may permit a more rigorous dis- 
cussion of finite-frequency noise with respect to pin- 
holes. Recently, first steps towards the discussion of 
timc-dcpcndcncc using this method have been made |42j . 
Also, electron-electron interactions describing the tradi- 
tional 1/f noise scenario for Josephson junctions should 
be included. 

We acknowledge stimulating discussion with John M. 
Martinis, W. Belzig and Yu. V. Nazarov. This work 
was financially supported by NSERC through a discov- 
ery grant and Quantum Works, EuroSQIP and Studiens- 
tiftung des deutschen Volkes. 
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